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Summary 

Severel energy- momentum "tensors" of gravitational field are considered 
and compared in the lowest approximation. Each of them together with 
energy-momentum tensor of point-like particles satisfies the conservation laws 
when equation of motion of particles are the same as in general relativity. It 
is shown that in Newtonian approximation the considered tensors differ one 
from the other in the way their energy density is distributed between energy 
density of interection (nonzero only at locations of particles) and energy 
density of gravitational field. 

Starting from Lorentz invariance the Lagrangians for spin-2, mass-0 field 
are constracted. They differ only by divergences. From these Lagrangians 
by Belinfante-Rosenfeld procedure the energy-momentum tensors are build. 
Only one of them is suitable for explaining the perihelion shift. This tensor 
does not coincide with Weinberg's one (directly obtainable from Einstein 
equation). 

It is noted that phenomenological field-theoretical approach (utilizing 
only vertices and propagators) can lead to modification of theory in the 
region of strong field, where till now no observational data are available. 

1 Introduction 



General relativity is a complete, elegant, and self-consistent theory. Yet there 
is a necessity to obtain gravity by field-theoretical means starting from flat 
spacetime, see e.g. [1-5]. It is widely believed that on this way even dropping 
the general covariance requirement we naturally get general relativity. It is 
supposed that in the lowest nonlinear approximation this is demonstrated in 
detail by Thirring [2]. Yet this conclusion can not be drawn from [2], see 
Sec.2. 

In this connection and also because the gravitational collapse is consid- 
ered as the greatest crisis in physics [6] the research into possible alternative 
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theories acquire especial significance. It is quite natural to make the first 
step and to consider the simplest processes by utilizing vertices; the graviton 
propagator is known by analogy with electrodynamics. 

In the lowest nonlinear approximation it is necessary to know only 3- 
graviton vertex. We assume the simplest possibility: the source of graviton 
is the energy-momentum tensor of two other gravitons. In higher approxi- 
mations probably other vertices will be needed. Along this path one can find 
out what theories are possible without assuming general covariance and a 
priori restriction on vertices. An important step in this direction was made 
by Thirring [1-2]. We continue his investigation in the same approximation 
and restrict ourselves to point-like classical particels as sources of gravitation. 
Mainly we are interested in the simplest system consisting of a Newtonian 
center and test particle moving in its field. 

In general relativity classical particles move along geodesies in Rieman- 
nian space. This is the incarnation of equivalence principle and it is more 
reliable than specific equation determining the gravitational field [9]. Be- 
sides, the equation of motion for particles are obtainable from variational 
principle without any appeal to general relativity, see equation (2) on page 
181 in [6]. It is easy to verify also that eq.(ll) in [2], derived in this way 
coincide exactly with that of general relativity. As to the equation deter- 
mining gravitational field, it is possible to think that the phenomenological 
field-theoretical approach will lead to more complicated algorithm for get- 
ting the field. An interesting possibility in this direction was pointed out by 
Schwinger [10], see also Sec. 6. 

It is reasonable to believe together with Einstein that for some reason 
or other the singular behaviour near the gravitational radius does not cor- 
respond to reality, see §15 in .Pais's book [11] and Einstein paper [12]. . 
At present the Schwarzschild singularity is concidered as fictitious by many 
researchers because the geometry is nonsingular there. See however the text 
after (2.2.6) in [13] and after eq.(9.40) in [14], where they say convincingly 
about physical singularity. By field-theoretical approach it is difficult to un- 
derstand why in a finite system the acceleration of a freely falling particle 
becomes unlimited when it nears the horizon. Such a behaviour should be 
connected with the fact that according to [15] the gravitational energy in 
vacuum outside the sphere of radius R goes to — oo for R — > r g . In confor- 
mity with this the energy of matter and gravitational field inside the sphere 
of radius R goes to +oo, in such away that total energy of spherical body 
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is equal to its mass. But if a theory predict that the absolute value of field 
energy outside sphere of radius R might be greater than total energy of a 
body then the analogy with electrodynamics suggests that the concept of 
external field becomes inapplicable [16]. The belief in general relativity in 
similar circumstances is based upon the concept of nonlocalazability of grav- 
itational energy, see, e.g. §20.4 in [6]. What is more, general relativity does 
not need as a rule the gravitational energy -momentum pseudotensor. 

The situation changes drastically, when we begin to constract gravity the- 
ory starting from flat spacetime. Here the gravitational energy-momentum 
tensor appears to be necessary to describe the 3-graviton interaction. The 
nonlinear correction to the motion of a test particle depends on the chosen 
energy-momentum tensor. The latter is build from field Lagrangian, which 
is not unique as one can add to it some divergence terms. This leads to dif- 
ferent energy-momentum tensors. They can give rise to gravitational energy 
densities, which may have even different signs. The question of sign of energy 
density is of interest by itself. Provided the sign turns out to be positive, one 
should expect the weakening of gravitational interaction at r ~ r g = 2GM 
in comparison with Newtonian one in order that the gravitational energy 
outside the sphere of radius r were much less than the mass M of the center. 
The possibility of decreasing of interaction at small distances is suggested 
also by the behaviour of attraction force between two bodies supported by 
Weil's strut, see, §35 [17]. 

In order to understand in what way the various energy-momentum tensors 
differ one from the other we consider the following tensors: Thirring's [1,2], 
Landau-Lifshitz's [16], Papapetrou- Weinberg's [9] and tensor obtained from 
the Lagrangian given in Exercise 7.3 in [6]. The second and third tensors 
are representatives of general relativity, the rest are build from Lagrangians 
of free field of spin-2, mass-0 particles and symmetrized by the Belinfante- 
Rosenfeld procedure, see, §1 Ch.7 in [18]. 

In the considered approach it is suitable to subdivide the 3-graviton vertex 
into three vertices in accordance with three possibilities for choosing two 
gravitons out of three to form the energy-momentum tensor,- the source for 
the third graviton, see Sec. 3. So three diagrams contribute to nonlinear 
correction to the field. In contrast with this the energy-momentum tensor, 
figuring in the solution of Einstein equation by iteration procedure, is so 
defined that the correction to field is given simply by means of propagator, 
i.e. by one diagram only, see Sec. 2. 
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The main result of the paper is this: starting from quadratic Lagrangians 
(differing by divergence terms) of spin-2, mass-0 particles, the energy- momentum 
tensors are constructed by Belinfante-Rosenfeld procedure. It turns out that 
only certain combination of these tensors is fitted for correct description 
of perihelion shift. This combination does not coincide with Papapetrou- 
Weinberg tensor. 

The investigation of possibilities of phenomenological approach to gravi- 
tation without use of general covariance seems to us very promising. Valuable 
undertaking in this direction was made in [19]. 

2 Thirring's energy- momentum tensor 

Throughout the paper we use 

9/u, = Vuv + h^ u (1) 

In Sections 2, 3, 5 rj^ = diag(l, — 1, — 1, — 1, ); in Sections 4 and 6 rj^ = 
diag(—l, 1, 1, 1, ). In this Sec. we use Thirring's notation [2]; both greek and 
latin indices run from to 3. 

The gravitational field is described by the symmetric tensor h^, which 
contains spin-2 and lower spins, see, e.g. [3]. The unnecessary spins (spin-1 
and one of spin-0) are excluded by Hilbert gauge : 

w\ v = (hr - \^ v h), v = o, h = h\, K = ^h. (2) 

One way to obtain Thirring's tensor is to start from general relativity La- 
grangian ^J—gR. If we remove terms with second derivatives of g^ v into 
divergence terms and drop the latter, we get the function G(x) in eq.(93.1) 
in [16]. Retaining in it only quadratic in h^ u terms we get 

G(x) = i[VA^' A - 2V> A "'" + 2h u ,V - h x h> x ] (3) 
This is equivalent to Thirring's Lagrangian [2] 

L= i[Vv,A^ ,A - 2^^^ + ~ txt X }- (4) 



4 



Here 

= - V/2/, V> = Vv", f = 87rG, G = 6.67 • lO^cm 3 ^ • sec 2 ' (5) 

Using instead of h^ u is justified because then the analogy with elec- 
trodynamics becomes more close: ip^ is analogous to vector- potential 
and has the same dimensionality, M\fG has the dimensionality of electro- 
magnetic charge. We note that the Lagrangian (4) exactly corresponds to 
Schwinger's Lagrangian [10], who uses the notation r]^ v = diag(—l, 1, 1, 1, ) 
and 2h% h = -hl u = -V- 

The canonical energy-momentum tensor following from (4) , has the form 

L=\ [ yv,AV^' A - \p,W X - 2</> V'l , (6) 

Using </?^ instead of ip^ u is handy as many expressions become more compact 
and the consequences of imposition of Hilbert gauge more clear. 

The energy-momentum tensor for a static point-like mass (Newtonian 
center) 

M 

T M „= M5(x)5, Q 5 v0 . (8) 
In linear approximation this source generate the field 

fM 

</V = - W 2 / = 4^]f| ^°^ ' h ^ u = 4< Mm°<W (9) 

satisfying Hilbert condition (2). For one Newtonian center = —GM/r. For 
several centers 

= -G£-^-. (10) 
n \r- r„l 



In terms of 



we have 



V/ = h ^ ~ h = h/, (11) 



V = 2 0<V; h = hf = -40 = -h. (12) 
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The energy density of field (9) is positive 

T°° = ^(V0) 2 . (13) 

/ 

But T 1& ought to be supplemented to symmetric one by the spin part: 

Q1& = ^ 7<5 + f 7<5. ( 14 ) 

For Newtonian center Thirring obtains 

t iS = "^(W) 2 M*o ( 15 ) 

So in this case 9 00 is negative 

r = - Jg (V0)2 - (16) 

Turning now to conservation laws of total energy-momentum we remind 

first how matters stand in general relativity. There the energy-momentum 

v 

tensor of point-like particles T^ u is connected with its counterpart in special 
relativity by the relation, see (33.4), (33.5) and (106.4) in [16]: 

v ijs 
^ T ^=T^ = Y,m a uV —5{x -x a {t)), u» = dx»/ds, (17) 

g is determinant of g^ v . In terms of the conservation laws are (see (96.1) 
in [16]) 

T% = [T^{r, TU + h TV )],, = \h^ u T^. (18) 

We shall see below that T^ T h TU can be interpreted as (part of) interaction 
energy-momentum tensor. 

As is known the equation of motion of particles in general relativity 

^ + r^"v=o (19) 

is contained in consevation laws. Indeed from 

T jk -j =f ik j + T j mj f mk + T k mj f jm = (20) 



p 

taking into account that from definition of T 3 in (17) 

t j \j = -\{-9)-H-9hT* + (-g)-iT*j, V mj = ±g, m , (21) 
we get 

T* d + T k mj T m3 = 0. (22) 
This is equivalent to (19), because [9] 

T> m , 3 = ^^5(x-x a (t)). (23) 

Going back to field-theoretical approach, we rewrite the equation of motion 
of particles (19) in the lowest approximation 

^ = ^ = ^t*V - h\,u^. (24) 

Just at such movement of particles the divergence of total energy-momentum 
ought to be zero, and inversely, from zero divergence follows eq. (24). From 
(23) and (24) we find 

T*„ = \h a /T^ - h 5 aa T^ , (25) 

This agrees with (22) in considered appoximation. With the same accuracy 
this can be rewritten as 

{r' 6 + r^h a 5 ) n = l -h^ 6 r^. (26) 

Using linearized Einstein equation for ip^ v 

f, v = T, v -^ u T , T = Tf , (27) 

we get 

T 75 , 7 = f^ 5 T aP = -h a /T^ = -lhafi*T#. (28) 
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Spin part of energy-momentum tensor is conserved by itself and do not con- 
tribute to conservation laws : 

T 75 , 7 = 0. (29) 

So from (26) and (28) it follows that the conserved tensor contains in 
itself the interaction tensor [2] 

T^ = T ia hJ. (30) 
But it is not symmetric. In order to understand the reason we have to 

S - 

consider the properties of T 7 in detail despite the fact that it does not take 
part in conservation laws written in the form of eqs.(26) and (28). According 
to known rules [2,18] we have 

f jk = _F iik ti - F kij ti - F ik \i , (31) 

f)r f 

F i*= (<p k a ri i f ) -<p i a ri k f ) ), L —L ■ (32) 

The antisymmetric part of (31) is contained only in the last term. For it we 
have 

-^%=f^-) V^-(/M ^ + (33) 

The first two terms on the right hand side symmetrize the interaction tensor, 
the last two terms symmetrize the canonical one. 

It is not seen directly from (6) and (31) that field energy-momentum 
tensor # 7<5 in (14) is symmetric. This agrees with the fact that the proof of 
symmetry utilizes the Euler-Lagrange equations for field which is considered 
as free [18]. We are interested in interacting field. So using linearized Einstein 
equation (27) with source, we get 

f(T a ^ k a - r a y a ) 

l(T ak K a -T a ih k a ). (34) 



dL 



<P a 



dL 



ak.i 



/(T^ fc a -T a V a ) 
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In the last equation we use the connection between ifj uu and h^, see (5). 
Substituting in (30) 7 — > j,5 — > /c, we see that the sum of (30) and (34) is 
symmetric. This result retains if we start from another Lagrangian differing 
from Thirring's one in (6) by divergence because the linearized equation 
remains the same. 

s 

One should take into account however that the symmetric part of T 3 can 
also contain terms of interaction type. So for the Lagrangian in (6) similarly 
to (34) we find 

_ F nk. _ F ki3 ti = [fT _ 2(/ ^v j + /[r a V« + ^VJ 

-2(^'*'V a + ^'*'V«,i) + 2^ , 'V„, (35) 

s ., 

As a result we get for T 3 

f 3k = 2[(<p i *> a i <p k a + ^ a ^ 3 a ) - <p ia , ia <p kj + ^' Q 'V fc a,< - V^'Va 

V fc, Va,i + ^V,,] - 2^ Q 'V<M + 2fT ja <p k a . (36) 
/ 

Here last but one term, added to T jk , makes it symmetric. The last term 
can be rewritten in terms of h^ v in the form, see (9) and (11), (30), 

-T 3a (h a k - ^r] a k h) = - f 3k + \r 3k h. (37) 

int T u 

So the symmetrization of T 3 in (30) is reduced to its replacement by ^T 3k h. 
This tensor is nonzero only where particles are present. For Newtonian cen- 
ters the corresponding energy density 

-T m h = -2T OO (38) 
2 

int 

is positive (contrary to our intuition and) contrary to T in (30), see (12) 
and (10), where h and <fi are given for Newtonian centers. 

We note that the use of linearized Einstein eq. (27) in the expression for 

s 

T 3 leads to that eq. (29) is satisfied only with considered accuracy. The 
presence of interaction energy-momentum tensor means the appearence of 
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such vertex: the energy-momentum tensor of matter together with one of 
gravitons serves as a source for other graviton, see Fig.l 

Thirring assumes that his tensor 9 1& (see (14), (6), (31)) is an analog of 
energy-momentum tensor figuring in the r.h.s. of Einstein equation when 
iteration procedure is used. In other words the nonlinear correction to field 
is given only by diagram (2a) in Fig. 2. On this Fig. the short straight line 
has only conditional meaning: it represents the source of gravitons, namely 
the energy-momentum tensor build of two gravitons (real or virtual) shown 
on Fig. 2 as joining the ends of this line. The graviton emerging from the 
middle of straight line is emitted or absorbed by this source. On diagram 
(2a) the energy-momentum tensor is build from gravitons of Newtonian cen- 
ter. On diagrams (2b) and (2c) one of the virtual gravitons of Newtonian 
center interact with energy-momentum tensor of two other gravitons. All 
three diagrams on Fig.2 correspond to one Feynman diagram obtained by 
contracting the short straight line to a point. 

The contribution to nonlinear correction for field from diagram (2a) is 
easy to obtain. Indeed, from (14), (6) and (15) we have 



Using now the field equation in Hilbert gauge with 9^ u from (16) and (39) 




(39) 



□ 




2 



(40) 



we find 



h 00 = -7(j> 2 , h lk = 



G¥ 



= 



GM 



r 



i, k = 1,2, 3. (41) 



Here easily verifiable relations 




4x l x k - 1 2 

Y7^ — 



(42) 



were used. The obtained h^ u satisfies Hilbert condition. Going over to h^ v = 
h^ v — \f]^vhi we find the following nonlinear corrections 




(43) 
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Index 2 in h indicating the order of correction in powers of G, is dropped 
for brevity. 

Finally from (12) and (43) we have 



ds 2 = g 00 dt 2 - (1 - 20 + 3(j) 2 )5 ik dx i dx k - : 



2 x %x fedoc doc 



(44) 



000 = 1 + 20 - 40 2 . 
The transition to spherical coordinates is given by the relations 



(45) 



5 ik dx i dx h = dr 2 + r 2 (d9 2 + sin 2 (fd(f 2 ), 



XiXkdx z dx k 



(dr) 2 . 



The nonlinear correction — 40 2 in goo m (45) is of special interest for us. The 
correct value necessary to explain the perihelion shift is +20 2 (if used in 
equation for geodesic, see. [19] and Sec. 6). The quantity g 00 is observable 
in principle by redshift [16] uj = c^o^oo)" 1 . Why the nonlinear correction 
in (45) turns out to be negative? It will appear later on that it is caused 
solely by the source (15), see eqs. (95) and (117). The sources (8) and 
(15) have different signs, but the corresponding fields have the same sign. 
The answer is simple. The correction in (45) is only a small (and negative) 
part. The larger and positive part goes for converting initially bare mass in 
Newtonian potential into a dressed one, see eq. (113). Now the negative sign 
of correction is clear: the mass of Newtonian center at infinitly large distance 
appears as M, but at finite distance the test particle feels a greater mass and 
greater attraction, because (15) is negative. 

There are some technical problems on how to take into account the con- 
tributions from diagrams (2b), (2c) and Fig.l. To simplify the problem as 
much as possible we may look only into corrections to Newtonian poten- 
tial. As it will appear in Sec. 6, it seems impossible to get the right value on 
any assumption about the contributions from these diagrams for the case of 
Thirring's tensor. 

We note here that Thirring obtained from his tensor the necessary cor- 
rection. Yet his result is objectionable as he used illdefined gauge 



□ 2 A = 



G 2 M 2 
4r 2 
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see eq. (83) in [2]. Namely the source of A fall of too slowly for large r and 
the integral in (112), see below, diverge for large r'. 

In the next two Sections we shall see how energy- momentum "tensors" of 
general relativity differ from Thirring's tensor. 



3 Landau-Lifshitz pseudotensor of energy- 
momentum 

This pseudotensor in the sense and approximation considered here is a tensor. 
In the lowest approximation with help of relation 

^gg ik « (1 + ^){jf k ~ h lk ) w rf k - h tk . (46) 
we get from eq. (96.9) in [16] 

fik _ 1 rTikJT ,m _ ZilZkm _ Tkm,pT ,i _ Tim,pT ,k , Tim,pTk , 
lQllG ' b ,l ,b : m ,L "mp 11 lb mp ~ 11 lb m,p 

\h™>% q > k - ^¥h' k + V lk (h mn , p h^ n + l -k m l m - -Ji pq , m h™> m )].{Al) 

Comparison with canonical tensor (6) shows that it is connected with f k by 
the relation 

t ik = > ik + ^ - h , k = _ 2f(fi ^ 
F tk = -^—;[h lkr h ln ' n - h u tl h kn , n - h kn *h n / + h in *h k n J. (48) 

From (14) we see that now in place of T %k stands F' lk . But T %k was a 
conserved quantity, see. (29). So F lk should rather play the role of interaction 
energy-momentum tensor. Indeed, taking into account that in considered 
approximation h ik satisfies the linearized Einstein equation 

h .3 _ h . 3 — h. 3 _i_ nrt U ,qr _ -IfW/^T (AO) 

we find 
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Now we check that coservation laws [16] 




are fulfilled. From (48), (28) and (50) we have 

t tk , = -\h^ k % q + h kn *T m - h>% k . (52) 
For matter energy-momentum tensor from (17) we get 

(- g ) f tk = V=gT lk « (1 + \)T lk , -gnl + h. (53) 

From here 

((-<?) T lk ) = T ik ,i + \h^ k , (54) 

the terms of order h 2 being dropped. Now it follows from (25), (52) and 
(54) that (51) is fulfilled. As seen from (53) here too there is a tensor, which 
is nonzero only where particles are located. Surprisingly it coincides with 
Thirring's interaction tensor, see (37) and text below it. 

p int 

Althouht —gT +t tk differs from Thirring's T tk + T + , for New- 
tonian centers they coincide. Now we turn to Newtonian approximation. 

According to Problem 1 in §106 in [16] the energy density of gravitational 

/ 

field in this approximation is given by T 00 in (13), (10). But there is also en- 
ergy density of interaction /Mp, where \i is density of particles. Using Poisson 
equation V 2 = inGfi and ignoring the problems connected with point-like 
nature of particles, we can write (utilizing integration by parts) 

J ^dV = - J ^(V0)W, (55) 

The density in the integrand on the r.h.s. contains now not only the energy 
density of interaction, but also the proper energy density of particle's self- 
field. The density on the l.h.s. is nonzero only at particle locations, the 
density on the r.h.s. is nonzero where the field is nonzero. The integration 
by parts deprive us the possibility to retain the previous physical meaning 
of integrand. If nevertheless we do this, then adding to (13) the energy 
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density in the r.h.s. of (55) we get the effective gravitational energy density 
in Newtonian approximation 

-^wr. (56, 

To bring this in agreement with t 00 we ought, according to a foot-note in 

v 

[16], take into account the contribution from (— g) T 00 - Let us do it. For t 00 
we have 

t°° = "^(V0) 2 , (57) 

57TG 

where is the potential of Newtonian centers. Now 

(-9) t 00 = v^T 00 « T 00 (l + |) « T 00 - 2«M, (58) 

see (12). The sought for agreement will be reached only after we rewrite a la 
Thirring [2] T 00 in terms of observables. From (17) 

T oo = ^m a ^5(x-x a (t))- (59) 
ds 

In the presence of gravitational field 

ds 2 = g 00 dt 2 (l-v 2 ). (60) 
Here v 2 is physical velocity, see §88 in [16]. Hence 



dx° 1 1 hoo > 

—— r^j - 

ds ^9m{l-v 2 ) ~ VI 2 



w 2 < 1. (61) 



Thus after going over to the observable velocity we obtain the term 

-^-H = -4>n. (62) 

detached from T 00 . Equation (12) was used to get the r.h.s.. Together with 
corresponding term in (58) this leads to 

-3^0 ^^(V0) 2 , (63) 
14 



where arrow corresponds to going over in (55) from integrand on the l.h.s. 
to the integrand on r.h.s.. Now the sum of (57) and (63) gives the expected 



The consideration of Newtonian approximation makes the following point 
of view very enticing: The energy density of an isolated point-like particle 
should be positive; Hilbert gauge exclude the unnecessary spins and then pos- 
itivity seems quite natural, because the presence of virtual gravitons should 
not make the energy density negative. The attraction is described by inter- 
action energy density and so the latter must be negative. Neither Thirring 
tensor nor Landau-Lifshitz one satisfies this requirement. The MTW ten- 
sor does. Unfortunately I failed to fit this idea into existing approach to 
gravitation, see Sec. 6. 

4 Papapetrou- Weinberg energy- momentum 
tensor. 

Einstein equation can be recast in such a way that gravitational energy - 
momentum "tensor" can be easily identified in coordinate system that goes 
over to Lorentzian at large distances from gravitating bodies [9]. In the 
lowest approximation this tensor have the form, see eq. (7.6.14) in [9]) 



VfMU = diag(-l, 1, 1, 1), 9vv = r]^ + V> V = K» = ~ h = ~Kv- ( 65 ) 



the first and second approximation in powers of h^ v . The indices are raised 
and lowered with rj^ v . In terms of h^ v we get 



(56). 



(64) 



In this Section we use Weinberg notation: 





2 (^V^,o- h, 



\K , (66) 
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h' a (-h --h --h ) + -h" (h a + h a 

+V^l\h x »h M - hh x x + ±h\M - hji% (67) 



R {2) = h x »{h KU >\ - \h vX /) + h\, v Jf\ K - \h\,M 

+ \Ka,xh KX ^ - -Ti^Ji^ - \hh KX M + h a k°. (68) 



For Newtonian center from (9-12) we obtain 
Nonzero components of t^ K are 



h v» = ~ h lu = -4<Mm)<W = -h T uv = -2(f)5„ u , h = h w = h T = -40 = -h. 

(69) 



too = "^(V0) 2 = U k = ^(75^ - U,V). (70) 

In Hilbert gauge from equation 

V 2 V = -16ttGV (71) 

we find, cf. with (40), (42), 

//on = ^f - = 30 2 , /,,, = -7G 2 M 2 ^. ( 72) 



It is easy to check that (72) satisfies the Hilbert condition (2). In terms of 
we have 

^oo = -20 2 , h lk = G 2 M 2 (^ - 7^). (73) 

In the expressions (71-73) h^ v is nonlinear correction. 

On the other hand in harmonic coordinates the Schwarzschild solution 
has the form [9] 



1 + r 2 



-dr 2 = -^^-dt 2 + (1 - <p) 2 (dx) 2 + r^z^_ dx i dx k (74) 
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So in the considered approximation this gives 

(?oo = -(l + 20 + 20 2 ), (75) 

g ik = l-2<p5 ik + (f) 2 {5 tk + ^). (76) 

From (69) we have Jiqq = —20, from (73) = — 20 2 , and there is agreement 
with (75). As to the nonlinear correction for g ik , in (73) it differs from the 
one in (76) by a gauge. Really, subtracting from h ik in (73) the nonlinear 
part of (76), we find 

G 2 M 2 ( ^ _ 8^ = ^2^^ + Afe , ); k% = ^ (7?) 

i.e. a gauge. 

Going back to t 0Q in (70), we note that this density is negative and does 
not coincide with any density of other tensors. At the same time the equa- 
tion of motion of particles is contained in the conservation laws of total 
energy-momentum tensor. We shall check it in considered approximation. 
For gravitational part the calculations give 

t^,. = -h\^ c + \h a T»° - l -h"^T pa - /^ A TV- (78) 

The energy-momentum tensor for particles, figuring in conservation laws, has 
a rather complicated form by construction [9] 

t m« = v ^ri KT gr7a g T p f a(B ~T MK + h\T aK + h\T a » 

« T^ K - -hT» K + W a T aR + h K a T c " i . (79) 

From here with the considered accuracy 

r"% =f "% + h^T™ + + h\T a ^ K . (80) 

So 

( t m« + t ^) K = f ^ k _ l -h^T pa + h^ K T aK + h a T^. (81) 
Further from (53) we get 

f » T"% - \h K T^. (82) 
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Taking into account (25) we see that the r.h.s. of (81) is zero. 

Now we turn to Newtonian approximation. Terms of interaction tensor 
are contained in both r^ v (three last terms in the r.h.s. of (79)) and in t^ K . 
From (64), (66-68), using (49), which preserve its form in the notation of this 
Section, we find the following terms of interaction tensor contained in t^ K : 

1 -1-1- 
~2 h ^ T ~ V^(h pa T pa - -hT) - -h% K . 

From here and the first equation in (70) we get in Newtonian approximation 

too = -^(V0) 2 -60T oo . (83) 

07TU 

Here is the same as in (10). From (79) and (53) we get in this approximation 
r °° = (-g)- l 2T 00 + 2h° T 00 » T 00 (l - -h) - 2/iooT 00 = T 00 + 60T 00 . (84) 

Thus in Newtonian approximation the interaction terms in the sum of (83) 
and (84) are cancelled out. The agreement with Newtonian approximation 
(56) is achieved in the same way as for Landau-Lifshitz tensor: T 00 on the 
r.h.s. of (84) detaches term (62), which is equivalent (in accordance with 
(63)) to 4^j(V0) 2 . Together with the first term on the r.h.s. of (83) this 
gives (56). 

Weinberg shows in detail that his energy-momentun tensor has all re- 
quiered characteristics. But this tensor does not help us to find energy- 
momentum tensor of two gravitons as represented by straight line on dia- 
grams of Fig.2. By construction Weinberg's tensor gives the gravitational 
field only via diagram (2a). The field-theretical description tell us that test 
particle is not quite passive. It does not simply follows the command "move 
along geodesic" but itself takes part in the creation of field in which it moves, 
see Fig.l and Fig. (2b, 2c). From this viewpoint one can expect that e.g. pho- 
ton and graviton scatter differently on Newtonian center as only in the latter 
case all three diagrams of Fig.2 contribute. 
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5 Misner-Thorne- Wheeler energy- momentum 
tensor 

In this Section it is handy for us to use again Thirring's notation. Up to 
divergence terms the Lagrangian (6) may be rewritten as [6] 

C = \f,u,x^' X - \ V ,^ X - ¥V W (85) 

The corresponding canonical energy-momentun tensor 

f BC 

T ]k = 7: <p^ k - ri Jk £, (86) 

acquires the form 

f, jk = fjk _ ^jkf ) fjk = ^kp^j _ Ifptfpj -2 V ju > k p ua ' a , T — —T — 2C. 

(87) 

Spin part is given by (31-32) with substitution L — > C. We dwell on differ- 
ences from Thirring's tensor. In symmetric in jk tensor 

piik + pkij = ^aij _ (p cuT^ (p k a + ^ai,k _ ^aa ^ki^ ^ 

-2^%<^ + (2^%7^ - ^ - ^ Va + tp ka ,«V ij + ^ a ^ ik (88) 

there is no derivatives over x l . This means that in divergence F^ lk ^ + F kl \i 
there are no terms of interaction tensor. In antisymmetric in jk tensor 

F iki = {^-^^ifPa-if^^H^ (89) 

such terms are present. Hence, using linearized Einstein equation (27) in the 
expression for F tk ^i, we get 

_ F ikJ i = f(Tt°<p k a - r fe Va) + ^A^a k ~ <p k a j ). (90) 

Terms with / together with (30) give symmetric interaction tensor in accor- 

/ 

dance with (34). Other two terms on the r.h.s. of (90) supplement T jk to 
symmetric one, see (87). So we get 

jk = f, jk+ f jk = ^k v ^3 _ l _ v ,k v ,j _ ^ a .^k _ pka ^ij 

-</^V<m - </^'Va,* + (y^' ,fc * + v ak ' j iW* + W</? J ' fe + 

2^^ k \ t - 2^\^ k \i + + <p aj > k )<p aa >° 

-2rf k {^\^\), - rf k C + /(T'V« - T ka ^ a ). (91) 
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From (91) and (30), using the relation between </? M „ and h^ u in (9), and taking 
into account (34), we find 

int 1 — — 1— — — — — — 

-{h ai ' j h h a,i + h ai ' k h j a,i) + (h aj ' k i + h ak \)h\ + ■>ir. a ,lr> h + -lir. J," '., 
-2^'%/^ + (h ak ' J + h ahk )K^ - 2rf k {h a \Ji\)^ - ^ k C] 

+i(T fcQ V +r^h a k ), (92) 

It is easy to verify that total energy-momentum tensor consisting of (17) and 
(92) is conserved. 

For Newtonian centers from (92) and (9), (12) we have 

^oo + f oo = 1 (V0)2 + 2M = r oo (g4) 

07TG 

For this system the MTW Lagrangian (93) coincide with Thirring's one. 
The same is true for canonical energy-momentum tensors, see (6) and (86- 
87), but spin parts are different. It follows from (91) and (87) that for 
Newtonian centers MTW spin part contributes only to interaction tensor, 
while Thirring's spin part contributes also to pure field part, see (15). 

We note now that in Hilbert gauge for static case (for Newtonian centers) 
the linearized Einstein equation can be written as 

V 2 ^ = -V 2 /^ = -lGnGT^, (95) 

see (116) below. As seen from (87) for this system T o — 0, i.e. there is no 
contribution to h 00 from diagrams of Fig.2. 

Comparing MTW and Thirring tensors in Newtonian approximation we 
see that addition of divergence terms to Lagrangian leads to the change in 
subdivision of energy density between purely field part and interaction part. 



6 On Schwinger's explanation of perihelion 
shift 

Schwiger gave an elementary explanation of perihelion precession [10]. His 
method is of interest in many respects. The offered algorithm for calculating 
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gravitational field is more in line with field-theoretical method than finding 
field via differential equation in geometrical approach. In linear approxi- 
mation both methods naturally agree. In higher approximation Schwinger's 
gravitational field, defined as a coefficient at 5T^ U in variation of amplitude 

6W(T) = J d A x5T^h*:\ 8T»\ V = 0, V = 2 ^f, (96) 

seems does not coincide with definition of h pu in (1), because the nonlinear 
correction to is \4> 2 , but not 2 , see the text below eq.(108). In the lowest 
approximation 

W(T) = 4-kG J d 4 xd A x'T^(x)D pi/pa (x - x')T^(x'), D^ pa = P^ pa D + (x - x') 

Pauper -r)ixvr)po). (98) 

The amplitude W is related to the energy of system E by the expression 

W{T) = - J dtE{t). (99) 

From the total energy E Schwinger singles out the interaction energy E mt . 
Nonlinear correction to the interaction energy of slowly moving planet with 
the Sun is given then by 

(2) r (7 3 r 

E mt = -GM / ^f mnt , (100) 

where t 00mt is taken in Newtonian approximation 

tmmt = _GMm v 11^^ 
4tt \x\ \x-R\ 

So the result is 

S" = ^(x), M.x) = ~. (102) 

Eq.(lOO) follows from eq.(4.34) in Ch.2 in [10] if one assums that energy- 
momentum tensor of gravitational field is given by 1/8 of the expression (15), 
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in which <fi is that of (10), i.e. the potential of Sun-planet system. The Sun's 
matter energy-momentum tensor is represented by (8). It is clear that (100) 
is given by contribution from diagrams (2b) and (2c). 

As shown in Sec. 3 the energy density of gravitational field may be cosidered 

(2) . 

as positive. In that case the expression for E mt in (102) changes sign, but 
then there is also a contribution from interaction energy density 

^ = + S(£ _ £)]. (103) 

R 

Using this instead of t 00mt in (100) and dropping the contribution from 5(x) 
(self-interaction is included in renormalization of Sun's mass) we get 

R 2 

This together with (102), taken with reversed sign, restores Schwinger's re- 
sult. We see that in this way the explanation of perihelion shift is achieved 
without the concept of negative energy density of gravitational field. 

So the nonlinear correction to the potential is taken into account by sub- 
stitution 

0-0 e// = 0(1 + ^0). (105) 

Therefore the Newtonian attraction is decreased by this correction. Qualita- 
tively this result is quite understandable: the correction term describes the 
interaction of the Sun with negative interaction energy and this corresponds 
to repulsion. The effective potential (105) produces the acceleration 

_iV/ = _,< (1 + ,) = _^ (1 _^). (106, 
ar r z r 

In general relativity the expression for acceleration of particle at rest 
depends upon the chosen coordinates. In Schwarzschild coordinates we have 

^= V ^ = 2GM > ( 107 ) 



see e.g. eqs. (A 61), (A 63) in Appendix in [13]. In harmonic coordinates 

2 



r = R — the expression takes the form 



GM GM GM . , ifto , 
= — (1 - — + •••)» (108) 



^(i + ^A-(^) 2 r 
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and does not contradict (106). We note also that in this coordinates goo is 
given by (75), i.e. the Newtonian potential is substituted by 0(1 + 0). 
The difference with (105) is caused by the fact that Schwinger obtained the 
correction to for use in Lagrangian method, cf. §106 in [16] and [21], not 
for correcting g Q0 . The relation is simple. Let us denote by |a and f3 the 
coefficients at 2 in corrections to in g 00 an d in Schwinger method. Then 
from eq. below (106.15) in [16] we have for particle at rest 



ds r ; — (a 1 

Loc- = v /l + ^oo-l+0+( v 2-2 

h 00 = 1 + 20 + a 



where 



So/?=f-i 



Now we shall see what is in store for us if we use the diagram method. 
First we dwell on properties of graviton propagator (98). The polarization 



factor P^ upa satisfies the relations 

iTP^T^ = t^T^ - hr = t» v % v = t^T^. (109) 
The scalar factor D + (x) has the representation 

D +{ x ) = T-M^ 2 ) = , n \o o 1 • , ( 110 ) 
and possesses the property 

/ dTD+{S ~ *> T) = W? I-oo (x - x'f- r 2 + »e = 4^hw[ (U1) 
For spherically symmetric body we have to deal with integrals of the kind 

d A x'D + {x-x')p{r') = — / = p(r') = - / dr'r' 2 p(r')+ / dr'r'p(r'), 

47T J y/x' 2 + X 2 — 2x ■ X' T JO Jr 

(112) 
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By the way it is seen from here that the derivative of Newtonian potential 
over r is determined only by the mass inside sphere of radius r. Assuming 
in (112) that p = we get 

IjfAVyO + jT ^VM = c(i-^). (113) 

Hence the divergent part at small r' appears only in the term, which is 
absorbed by mass renormalization. 
So the source (13) generate the field 

h 00 = IQtiG J d 4 x'D + (x - x')T 00 (x') -(f) 2 . (114) 

The arrow shows that the divergent part is included in mass renormalization. 
The r.h.s. of (114) can be obtained also from the solution of wave equation 
derived from (114) by action of the operator d 2 = V 2 — and taking into 
account that 

-d 2 D + {x - x') = S(x - x') , V 2 «^. (115) 

We note also that 

V = 2h s * = lQnG j d A x'D + (x - x')% v (x'). (116) 

First we try MTW energy-momentum tensor as a source of gravitational 
field. The canonical part of MTW tensor in Hilbert gauge has the form 

T 7<5 = V' 7 - \h> 6 h>~< 

-l 7 ^^ - \h' X kx}} = T 1& - \^ 5 f. (117) 

For Newtonian centers this expression coincides with Thirring's one, see 
eq.(6). As seen from (117) in this static case T 00 = 0. So there is no 
contribution to /loo from this source. As noted in Sec. 5 the nonlocal part of 

s 

T ^ u is zero for Newtonian centers. So assuming, as in the case of Weinberg 
tensor, that the contribution from diagrams of Fig.l and (2b), (2c) should 
not be taken into account, we cannot obtain g o necessary for explaining the 
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perihelion shift; even a combination of MTW and Thirring tensors will not 
help in this case. 

Thus we have to take into account all diagrams. As we know from general 
relativity, where the vertices are known, the extracting post-Newtonian cor- 
rections from scattering amplitudes is a cumbersome and labour consuming 
job [21, 22]. For this reason we assume here the Schwinger method. 

Returning to Thirring tensor we find that the total correction to <fi from 
all diagrams of Fig. 1-2 is zero. Indeed starting from the energy density (56), 
Schwinger obtained the correction to as |0 2 . Therefore the corresponding 
contribution from spin part (15) is 4</> 2 . This correction comes from diagrams 
(2b) and (2c). Similarly to (103), (104) we find that the correction to <fi from 
diagram of Fig.l is — 20 2 . Finally the diagram (2a), treated according to 
Schwinger method, contributes — 2(f) 2 . So we get zero instead of |0 2 , necessary 
for use in Lagrangian method, see §106 in [16]. 

Now it easy to see that a linear combination of MTW and Thirring tensors 
with weights \ and | gives the necessary correction |0 2 . As the MTW tensor 
alone is unable to explain the perihelion shift, we have to give up the tempting 
disire to consider the gravitational energy density of an isolated point-like 
center as positive. 

Besides the nonlinear correction (102) there is a nonrelativistic linear cor- 
rection to the Sun-planet interaction energy. All nonrelativistic corrections 
have the form [10] 

§£v + ^V, V = v*. (118) 

Relativistic corrections may be treated either according to Schwinger or by 
using (118) together with Newtonian potential energy V in relativistic equa- 
tion of motion in flat space, see i.g. Ch.5, §1 in [20]. Both methods give the 
correct result. An elaborate analysis of different aspects of planet motion 
in this approximation is given in [19]. Finally we remind that three other 
famous observational effects of general relativity are explained already by 
linear approximation [2,10,19]. 



7 Concluding remarks 

Though general covariance was not assumed, the gauge invarians is of course 
retained [2,10]. For this reason the weak gravitatinal waves in flat space are 
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described as in general relativity. All considered tensors give the same energy- 
momentum tensor for the plane gravitational wave. There are no a priori 
reasons to believe that field-theoretical approach will give the same result 
as general relativity. Similar assertion was made in [21]. Our investigation 
shows that there is still much to be done to synthesize the geometrical and 
field-theoretical aspects of gravitations. 

I wish to thank V.I.Ritus, S.L.Lebedev, I.V.Tyutin, and M.A.Vasiliev for 
useful discussions. I am particularly indebted to V.Ya. Fainberg for remarks 
contributing to improvement of the manuscript. This work was supported 
by the Russian Fund for Fundamental Research (Grants No 96-02-17314 and 
96-15-96463). 

Figure captions 

Fig.l: The second rank tensor formed from matter energy- momentum 
tensor and graviton is a source for another graviton. 

Fig. 2: 3-graviton vertex. Short straight line serves only to distinguish the 
roles of participating gravitons: energy-momentum tensor is formed from two 
gravitons joining the straight line at its ends, this energy-momentum tensor 
serves as a source for graviton emerging from the middle of the straight line. 
Crosses represent external field sources. 
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